Search for: All records

The Northern and Southern Patagonian Icefields are rapidly losing volume, with current volume loss rates greater than 20 km³a⁻¹. However, details of the spatial and temporal distribution of their volume loss remain uncertain. We evaluate the rate of 21st-century glacier volume loss using the hydrological balance of four glacierised Patagonian river basins. We isolate the streamflow contribution from changes in ice volume and evaluate whether the rate of volume loss has decreased, increased, or remained constant. Out of 11 glacierised sub-basins, seven exhibit significant increases in the rate of ice volume loss, with a 2006–2019 time integrated anomaly in the rate of glacier volume loss of 135 ± 50 km³. This anomaly in the rate of glacier-volume-loss is spatially heterogeneous, varying from a 7.06 ± 1.69 m a⁻¹increase in ice loss to a 3.18 ± 1.48 m a⁻¹decrease in ice loss. Greatest increases in the rate of ice loss are found in the early spring and late summer, suggesting a prolonging of the melt season. Our results highlight increasing, and in some cases accelerating, rates of volume loss of Patagonia's lake-terminating glaciers, with a 2006–2019 anomaly in the rate of glacier volume loss contributing an additional 0.027 ± 0.01 mm a⁻¹of global mean sea-level rise.

 

Tropical glacier melt provides valuable water to surrounding communities, but climate change is projected to cause the demise of many of these glaciers within the coming century. Understanding the future of tropical glaciers requires a detailed record of their thicknesses and volumes, which is currently lacking in the Northern Andes. We calculate present-day (2015–2021) ice-thicknesses for all glaciers in Colombia and Ecuador using six different methods, and combine these into multi-model ensemble mean ice thickness and volume maps. We compare our results against available field-based measurements, and show that current ice volumes in Ecuador and Colombia are 2.49 ± 0.25 km³and 1.68 ± 0.24 km³respectively. We detected no motion on any remaining ice in Venezuela. The overall ice volume in the region, 4.17 ± 0.35 km³, is half of the previous best estimate of 8.11 km³. These data can be used to better evaluate the status and distribution of water resources, as input for models of future glacier change, and to assess regional geohazards associated with ice-clad volcanoes.

 

Atmospheric and oceanic warming over the past century have driven rapid glacier thinning and retreat, destabilizing hillslopes and increasing the frequency of landslides. The impact of these landslides on glacier dynamics and resultant secondary landslide hazards are not fully understood. We investigated how a 262 ± 77 × 106 m3 landslide affected the flow of Amalia Glacier, Chilean Patagonia. Despite being one of the largest recorded landslides in a glaciated region, it emplaced little debris onto the glacier surface. Instead, it left a series of landslide-perpendicular ridges, landslide-parallel fractures, and an apron of ice debris—with blocks as much as 25 m across. Our observations suggest that a deep-seated failure of the mountainside impacted the glacier flank, propagating brittle deformation through the ice and emplacing the bulk of the rock mass below the glacier. The landslide triggered a brief downglacier acceleration of Amalia Glacier followed by a slowdown of as much as 60% of the pre-landslide speed and increased suspended-sediment concentrations in the fjord. These results highlight that landslides may induce widespread and long-lasting disruptions to glacier dynamics. 

Ancient glaciated landscapes reveal interactions among ice dynamics, bed composition, and preglacial river networks.

Subglacial landscapes, revealed in regions of recent ice-sheet retreat, provide a window into ice-sheet dynamics and interactions with evolving subglacial topography. Here, we document landscape evolution beneath the southern Laurentide Ice Sheet of North America since the end of the Pliocene, 2.6 million years (Ma) ago, by reconstructing the isostatically adjusted preglacial surface and modern bedrock topography at 250 m horizontal resolution. We use flow routing to reconstruct drainage networks and river longitudinal profiles, revealing the pattern and extent of their glacially forced reorganization. The overall mean Quaternary (2.6 Ma ago to present) erosion rate is 27 m/Ma, rising within ice-streaming corridors to 35 m/Ma (and locally reaching 400 m/Ma) and falling to 22 m/Ma in non–ice-streaming regions. Our results suggest that subglacial erosion was sufficient to lower the southern Laurentide Ice Sheet into warmer environments, thereby enhancing ablation and reducing ice-sheet extent over time. 

Abstract. Depressions – inwardly draining regions – are common to many landscapes. When there is sufficient moisture, depressions take the form of lakes and wetlands; otherwise, they may be dry. Hydrological flow models used in geomorphology, hydrology, planetary science, soil and water conservation, and other fields often eliminate depressions through filling or breaching; however, this can produce unrealistic results. Models that retain depressions, on the other hand, are often undesirably expensive to run. In previous work we began to address this by developing a depression hierarchy data structure to capture the full topographic complexity of depressions in a region. Here, we extend this work by presenting the Fill–Spill–Merge algorithm that utilizes our depression hierarchy data structure to rapidly process and distribute runoff. Runoff fills depressions, which then overflow and spill into their neighbors. If both a depression and its neighbor fill, they merge. We provide a detailed explanation of the algorithm and results from two sample study areas. In these case studies, the algorithm runs 90–2600 times faster (with a reduction in compute time of 2000–63 000 times) than the commonly used Jacobi iteration and produces a more accurate output. Complete, well-commented, open-source code with 97 % test coverage is available on GitHub and Zenodo. 

Abstract. We present Glacier Image Velocimetry (GIV), an open-source and easy-to-use software toolkit for rapidly calculating high-spatial-resolutionglacier velocity fields. Glacier ice velocity fields reveal flow dynamics, ice-flux changes, and (with additional data and modelling) icethickness. Obtaining glacier velocity measurements over wide areas with field techniques is labour intensive and often associated with safetyrisks. The recent increased availability of high-resolution, short-repeat-time optical imagery allows us to obtain ice displacement fields using“feature tracking” based on matching persistent irregularities on the ice surface between images and hence, surface velocity over time. GIV isfully parallelized and automatically detects, filters, and extracts velocities from large datasets of images. Through this coupled toolchain and aneasy-to-use GUI, GIV can rapidly analyse hundreds to thousands of image pairs on a laptop or desktop computer. We present four example applicationsof the GIV toolkit in which we complement a glaciology field campaign (Glaciar Perito Moreno, Argentina) and calculate the velocity fields of smallmid-latitude (Glacier d'Argentière, France) and tropical glaciers (Volcán Chimborazo, Ecuador), as well as very large glaciers (Vavilov Ice Cap,Russia). Fully commented MATLAB code and a stand-alone app for GIV are available from GitHub and Zenodo (see https://doi.org/10.5281/zenodo.4624831, Van Wyk de Vries, 2021a). 

Abstract Landslides pose a major natural hazard, and heterogeneous conditions and limited data availability in the field make it difficult to connect mapped landslide inventories to the underlying mass-failure mechanics. To test and build predictive links between landslide observations and mechanics, we monitored 67.89 h of physical experiments in which an incising and laterally migrating river generated landslides by undercutting banks of moist sand. Using overhead photos (every 20 s) and 1-mm-resolution laser topographic scans (every 15–30 min), we quantified the area, width, length, depth, volume, and time of every visible landslide, as well as the scarp angles for those within 3 min prior to a topographic scan. Both the landslide area–frequency distribution and area–volume relationship are consistent with those from field data. Cohesive strength controlled the peak in landslide area–frequency distribution. These results provide experimental support for inverting landslide inventories to recover the mechanical properties of hillslopes, which can then be used to improve hazard predictions. 

Abstract. Depressions – inwardly draining regions of digital elevation models – present difficulties for terrain analysis and hydrological modeling. Analogous “depressions” also arise in image processing and morphological segmentation, where they may represent noise, features of interest, or both. Here we provide a new data structure – the depression hierarchy – that captures the full topologic and topographic complexity of depressions in a region. We treat depressions as networks in a way that is analogous to surface-water flow paths, in which individual sub-depressions merge together to form meta-depressions in a process that continues until they begin to drain externally. This hierarchy can be used to selectively fill or breach depressions or to accelerate dynamic models of hydrological flow. Complete, well-commented, open-source code and correctness tests are available on GitHub and Zenodo. 

Climatic warming following the Last Glacial Maximum caused the southern Laurentide Ice Sheet (LIS) to begin ∼2,000‐year cycles of retreat and readvance whose cause remains ambiguous. By developing a marine‐calibrated chronology of southern LIS position, we counterintuitively demonstrate that between 17.6 and 11.3 ka, ice advanced during times of northern‐hemisphere warming and retreated during times of northern‐hemisphere cooling. Here we propose a cyclical feedback: Meltwater from ice retreat cooled the northern hemisphere by weakening the Atlantic Meridional Overturning Circulation (AMOC). This eventually lead to ice‐sheet readvance, which reduced and rerouted meltwater discharge, and thereby allowed the AMOC to strengthen and the northern hemisphere to warm. Our data suggest that this antiphased ice–climate interaction, paced by ice‐sheet response time, was initiated by synchronous warming and ice retreat ∼18.7–17.6 ka (corresponding to the Erie “Interstade”) and reached its apex during the Younger Dryas.

 

Abstract. Calculating flow routing across a landscape is a routine process in geomorphology, hydrology, planetary science, and soil and water conservation. Flow-routing calculations often require a preprocessing step to remove depressions from a DEM to create a “flow-routing surface” that can host a continuous, integrated drainage network. However, real landscapes contain natural depressions that trap water. These are an important part of the hydrologic system and should be represented in flow-routing surfaces. Historically, depressions (or “pits”) in DEMs have been viewed as data errors, but the rapid expansion of high-resolution, high-precision DEM coverage increases the likelihood that depressions are real-world features. To address this long-standing problem of emerging significance, we developed FlowFill, an algorithm that routes a prescribed amount of runoff across the surface in order to flood depressions if enough water is available. This mass-conserving approach typically floods smaller depressions and those in wet areas, integrating drainage across them, while permitting internal drainage and disruptions to hydrologic connectivity. We present results from two sample study areas to which we apply a range of uniform initial runoff depths and report the resulting filled and unfilled depressions, the drainage network structure, and the required compute time. For the reach- to watershed-scale examples that we ran, FlowFill compute times ranged from approximately 1 to 30 min, with compute times per cell of 0.0001 to 0.006 s.